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Abstract — The Hyrcanian forests are green belt stretching 
over the northern slopes of the Alborz mountain ranges and 
cover the southern coasts of the Caspian Sea. The climate of this 
region is controlled by several components of a regional 
atmospheric circulation pattern and is strongly modulated by a 
complex topography and the maritime effect of the Caspian Sea. 
Climate change will accelerate the hydrologic cycle, altering 
rainfall, and the magnitude and timing of runoff. Hyrcanian 
forests might become one of the most vulnerable areas in the 
world regarding climate change. Therefore, the purpose of this 
paper is to assess the impacts of climate change on surface runoff 
from the Hyrcanian forests in the North of Iran. To study the 
effects of climatic variations, the SWAT model was implemented 
to simulate the hydrological regime and the SUFI-2 algorithm 
was used for parameter optimization. The climate change 
scenarios were constructed using outcomes of three General 
Circulation Models (CGCM2, HadCM3, and SCIR02) for three 
emission scenarios (A1F1, A2 and Bl) by adjusting the baseline 
climatic variables that represent the current precipitation and 
temperature patterns. The study results for 2040-2069 
compared with the present climate showed changes in surface 
runoff by -1.3%, 5% and -1.2% for the A1F1, A2 and Bl 
scenarios, respectively. Monthly variations show pronounced 
increases in discharge in the wet season (February-May) and 
decrease in dry season (July-September). The results highlight 
the strong impact of climate change in surface runoff and reflect 
the importance of incorporating such analysis into adaptive 
management. 

Index Terms — Climate Change, Hyrcanian Forests, Surface 
Runoff, SWAT, SWAT -CUP, Iran. 

I. INTRODUCTION 

Hyrcanian forests stretch out from sea-level up to an altitude 
of 2,800 m and encompass different forest types by the virtue 
of their 80 different woody species (trees and shrubs). The 
area is rich in hardwood species, but there are only four 
genera of endemic softwood (conifer) trees including yew, 
Greek juniper, oriental arbor-vitae and Italian cypress. 
However, based on the studies of Fadaiey Khojasteh et al. 
(2010) three genera of Mesozoic Gymnosperms were 
recognized. The primary function of the Hyrcanian forests, 
other than wood production, is supportive and environmental. 
They play a vital role in the conservation of soil and water 
resources and keep nature at balance on these susceptible 
steep mountain slopes. However, rapid urbanization and 
industrialization, intensive grazing, over-utilization of forests 
for firewood production and farming is destabilizing the 
forest and the environemnts around it. 
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Over the last few decades, swift forest degradation has 
brought about a number of environmental, social and 
economic impacts including soil erosion, floods, degradation 
of farmlands and habitats, reduction of biodiversity and 
natural resources, and air and water pollution. Furtheremore, 
manipulation of forest ecosystems has threatened a number of 
animal species such as fallow deer, roe deer, wolf, fox, wild 
cat, leopard, pheasant and trout. 

In recent years, climate change is one of the most important 
phenomena that threatens this unique ecosystem. The 
consensus of atmospheric scientists is that the earth is 
warming, and as global temperatures increase, the hydrologic 
cycle is becoming more vigorous. The IPCC has reported that 
there has been a very likely increase (probability 90-99%) in 
precipitation during the 20th century in the mid-to-high 
latitudes of the Northern Hemisphere. According to the 
Fourth Assessment Report (AR4) of IPCC, global mean 
surface temperature, precipitation and extreme events such as 
heavy precipitation and droughts have changed significantly, 
and the changes are very likely to continue (IPCC 2007).The 
rises of earth near-surface air temperature and changes in 
precipitation patterns are prominent features of climate 
change; these two factors impact almost all other hydrological 
processes. All Atmospheric-Ocean General Circulation 
Models (AOGCMs) predict a rise in earth surface temperature 
and rainfall intensity and amount due to increasing in 
greenhouse gasses (GHG) concentration over the coming 
century (Kami et al. 2010). 

A warmer climate will accelerate the hydrologic cycle, 
altering rainfall, magnitude, and timing of runoff. Warm air 
holds more moisture and increases evaporation of surface 
moisture. With more moisture in the atmosphere, rainfall and 
snowfall events tend to be more intense, increasing the 
potential for floods (Dhar and Mazumdar 2009). Using 
present day precipitation patterns, studies have shown that 
higher temperatures lead to increased evaporation rates, 
reductions in surface runoff, and increased the frequency of 
droughts (Ficklin et al. 2009). The changes in flow 
characteristics resulting from climate change depend on 
individual catchment characteristics. In particular, basin 
geology and elevation are first-order controls on the timing 
and magnitude of basin runoff to climate change (Hamlet and 
Lettenmaier 2007). Nearly all regions of the world are 
expected to experience a net negative impact of climate 
change on water resources. But the intensity and 
characteristics of the impact, however, can vary significantly 
from region to region (Abbaspour et al. 2009). Reliable 
predictions of the quantity and rate of runoff are needed to 
help decision makers in developing watershed management 
plans for better soil and water conservation measures. 

Many recent studies have focused on the potential effects of 
climate change on water resources including water quality and 
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quantity. Gosain et al (2006) simulated the impacts of a 
2041-2060 climate change scenario on stream discharges 
from 12 major river basins in India, ranging in size from 1,668 
to 87,180 km2. Stream discharge was found to generally 
decrease, and the severity of both floods and droughts 
increased in response to the climate change projection. Aimed 
to the prediction of surface runoff in the upper Mississippi 
River basin, Jha et al (2006) used various global climate 
models to predict surface runoff in the upper Mississippi 
River basin. Study results showed a wide range of changes, 
from a 6% decrease to a 51 % increase depending primarily on 
precipitation patterns. Abbaspour et al (2009) used the 
hydrologic program Soil and Water Assessment Tool 
(SWAT) (Arnold et al. 1998) to study the impact of future 
climate on water resources availability in Iran. Future climate 
scenarios for periods of 2010-2040 and 2070-2100 were 
generated from the Canadian Global Coupled Model (CGCM 
3.1) for scenarios A1B, Bl, and A2. Analysis of daily rainfall 
intensities indicated more frequent and larger-intensity floods 
in the wet regions and more prolonged droughts in the dry 
regions. Chang and Jung (2010) estimated potential changes 
in annual, seasonal, and high and low runoff and associated 
uncertainty in the 218 sub-basins of the Willamette River 
basin of Oregon. The seasonal variability of runoff is 
projected to increase consistently with increases in winter 
flow and decreases in summer flow. Zarghami et al (2011) 
used LARS-WG and General Circulation Models (GCM) 
outputs for prediction the climate change on the East 
Azerbaijan Province in Iran. The research outcomes using the 
artificial neural network (ANN), showed dramatic reductions 
in the flows. Azari et al (2015) simulated the impacts of a 
2040-2069 climate change impacts on surface runoff in 
Gorganroud river basin in the North of Iran. The study results 
showed an increase in annual surface runoff of 5.8%, 2.8% 
and 9.5% and an increase in sediment yield of 47.7%, 44.5% 
and 35.9% for the A1F1, A2, and Bl emission scenarios, 
respectively. 

The above studies indicate that watershed processes are very 
sensitive to changes in precipitation and temperature and can 
vary significantly from region to region. Therefore, 
quantifying Hydrological impacts of Climate Change and 
future conditions will be valuable in understanding and 
predicting discharge processes as well as watershed-scale 
sustainable water management. The potential future changes 
in sediment load Also should be seen as an important 
requirement for sound river basin management. In this study, 
we evaluate the potential impacts of climate changes on 
surface runoff in the Hyrcanian forests in the north of Iran. 
For this purpose, we used SWAT to simulatethe surface 
runoff with three AOGCS climate models (CGCM2, CSIRO, 
and HadCM3) for the time period of 2040-2069 under A1F1, 
A2 and Bl greenhouse gas emissions scenarios. This paper 
contributes to the scientific understanding of changing surface 
runoff in this region. 

II. Material and methods 

2.1 Study Area 

The Hyrcanian forest stretches from Astara in the northwest to 
the Gorgan vicinity in the northeast of Iran. This area is 
approximately 800 km long and 110 km wide and has a total 


area of 18,500 km 2 comprising 15 % of the total Iranian 
forests and 1.1 % of the country’s area. 

The study region located between 35° 47-36° 35' N and 50° 
34-54° 10' E in the Mazandaran province (Figure 1). 
Agriculture, forests, and range lands dominants the land use. 
The elevation ranges from -26m at the outlet to 5595 m at the 
top of Damavand peak in the south of the area. The annual 
rainfall varies from 231 mm to 1200 mm. The minimum and 
the maximum temperatures in the province ranges from 9°C 
and 18.1°C, respectively. The temperature of the warmest 
month ranges from 28 to 35 °C while that of the coldest month 
is between 1.5 and 4 °C. Summer temperature ranges between 
20 and 30 °C. 



Fig. 1 Location of the part of Hyrcanian forests in the 
Mazandaran province 


2.2 The SWAT Model 

The Soil and Water Assessment Tool (SWAT) (Arnold et al. 
1998) is a physical process based model to simulate 
continuous-time landscape processes at a catchment scale. In 
SWAT watershed is divided into hydrological response units 
(HRUs) based on soil type, land use and slope classes that 
allow a high level of spatial detail simulation. The major 
model components include hydrology, weather, soil erosion, 
nutrients, soil temperature, crop growth, pesticides 
agricultural management and stream routing. The model 
predicts the hydrology at each HRU using the water balance 
equation, which includes daily precipitation, runoff, 
evapotranspiration, percolation and return flow components. 
The surface runoff is estimated in the model using two options 
(i) the Natural Resources Conservation Service Curve 
Number (CN) method and (ii) the Green and Ampt method. 
The percolation through each soil layer is predicted using 
storage routing techniques combined with the crack-flow 
model. The evapotranspiration is estimated in SWAT using 
three options (i) Priestley-Taylor, (ii) Penman-Monteith and 
(iii) Hargreaves. The flow routing in the river channels is 
computed using the variable storage coefficient method, or 
Muskingum method (Arnold et al. 1998). The wide range of 
SWAT applications underscores that the model is a very 
flexible and robust tool that can be used to simulate a variety 
of watershed problems. Hence SWAT was selected for use in 
this study because of its ability to simulate regional water flow 
at a watershed scale and to provide effective results. 
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2.3 Data and Model setup 


2.5 Future Climate Data 


Land use map extracted from the interpretation of Land sat 
TM (30m resolution) satellite imagery, based on field 
investigation, which contains seven different land use 
classes. Soil map and texture was obtained from Iranian 
ministry of Agriculture which has a spatial resolution of 
1:250,000 and includes a set of estimated physical and 
chemical soil properties. The catchment area of the 
Mazandaran province was delineated and discretized into 
sub-basins using a 30m DEM. Daily observed climate data 
including daily precipitation and temperature were obtained 
for 25 stations from the Iranian Meteorological Organization 
and the Water Resources Management Organization 
(WRMO) of Iran. Daily river discharge data required for 
calibration-validation were obtained from the WRMO of 
Iran. The monthly discharge data from 20 hydrometric 
stations within the basin for 34 years were used for model 
calibration and validation. Three slope classes including 0- 
15, 15-30, 30-60 were used in HRU definition. With these 
specifications, a total of 372 sub basins and 2535 HRUs 
were delineated in the study area. 


2.4 Calibration and sensitivity analysis 


Parameter optimization and uncertainty analysis were done 
using the Sequential Uncertainty Fitting Program SUFI-2 
(Abbaspour, 2007). In this algorithm, all uncertainties 
(parameter, conceptual model, input, etc.) are mapped onto 
the parameter ranges as the procedure tries to capture most of 
the measured data within the 95% prediction uncertainty 
(95PPU). Two indices were used to quantify the goodness of 
calibration/uncertainty performance. Two indexes define the 
strength of calibration and the prediction uncertainty: 
P-factor, which is the percentage of data bracketed by the 
95PPU band (maximum value 100%), and the R-factor, which 
is the average width of the band 95PPU divided by the 
standard deviation of the corresponding measured variable. 
Model evaluation is an essential measure to verify the 
robustness of the model. The performance of the model for 
simulating discharge is evaluated by Nash-Sutcliffe efficiency 
(E ns ) (Eq. 1), and the coefficient of determination (R 2 ) (Eq. 
2). E ns ranges from negative infinity to 1, with 1 denoting a 
perfect model agreement with observation (Nash and Sutcliffe 
1970). 
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In these equations, n is number of observed data, ^ and 


are observed and simulated data, respectively, on each 
time step i (e.g., day or month), and 3^ si7n are mean 

values for observed and simulated data, respectively. We 
considered 1972-1996 and 1997-2006 as the simulation 
periods for calibration and validation, respectively. The first 
two years was considered as a warm-up period in which the 
model was allowed to initialize and approach reasonable 
initial values for model state variables. 


A common approach for assessing future runoff conditions is 
to use climate model projections in combination with 
hydrological models. In this study, we used data from climate 
simulations statistically downscaled by the Climatic Research 
Unit, University of East Anglia. The three Global Climate 
Models (GCM) used were: CGCM2 (Coupled Global Climate 
Model) from Canadian Center for Climate Modeling and 
Analysis, HadCM3 from Hadley Centre for Climate 
Prediction and Research and SCIR02 from Australia's 
Commonwealth Scientific and Industrial Research 
Organization. Scenarios with the highest (A1FI scenario - 
970 ppm by 2100), lowest (B1 scenario - 550 ppm by 2100) 
and plausible (A2 scenario - 845 ppm by 2100) projected 
C02 concentrations were chosen for this study. Monthly 
maximum temperature, minimum temperature, and 
precipitation on a 0.5° grid are available for globe from 2001 
to 2100 (Mitchell et al. 2004). 

Climate change scenarios were developed using downscaled 
monthly average total precipitation and monthly mean 
temperature data. The baseline data was from 1971-2000. 
Initially, the GCM gridded data were spatially interpolated to 
the target stations using inverse distance weighted averaging 
of four native neighbors. Taking the center as the grid point 
for each grid box, we used. 



Where Si is the downscaled site-specific GCM projection at 
site i, p k is the GCM projection at the cell k , d i>k is the distance 
between site i and the center of cell k , m -3 is used in this study 
(Liu and Zuo, 2012). Then Change Factor (CF) method was 
used to generate climate change scenarios for 2040-2069. The 
CF method involves adjusting the observed daily temperature 
(T 0 bs,d) by adding the difference in monthly temperature 
predicted by the climate model (GCM or RCM) between the 
future and the reference period (T CM5 f ut5m T CM?ref?m ). To obtain 
daily temperature at the future horizon (T ad j ?fut?d ) we used Eq. 
4. The adjusted daily precipitation for the future horizon 
(Padj,fut,d) is obtained by multiplying the precipitation ratio 
(PCM,fut,n/P cm, ref,m) with the observed daily precipitation 
(Pobs’d) (Eq. 5)(Chen et al. 2011). 

^'adj f fut f d ^obs.,d ^ (.TcMJutym ^ CM f ref f m ) 

( 4 ) 

Pq dj r futfd Pobs.d ^ (.PcM.,fut,m PcM r rsf r 7n) 

( 5 ) 

Finally, daily data for future climate projections by GCMs 
under different greenhouse gas emissions scenarios for every 
station were used as inputs to the modified SWAT to project 
the watershed-scale changes in hydrological components in 
the 2040-2069. 


III. Results and discussion 

3.1 Model calibration and verification 

The SWAT model was calibrated based on daily measured 
discharge at 20 stations within the watershed. First, 
Sensitivity analysis Using SUFI-2 in SWAT-CUP was 
performed to evaluate the effect of parameters on the 
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performance of SWAT in the simulating runoff. So sensitivity 
analysis, Calibration, and validation of SWAT model were 
done for every station separately. Figure 2 compares 
graphically measured and simulated monthly surface runoff 
with 95% prediction uncertainty band for the calibration and 
validation period at 8 stations located in main outlets. 


In addition to the visual comparison, the ststistics of the 
results for the eight discharge stations above are given in 
Table 1. The overall N S e and R2 for the calibration and 
validation periods indicated a close relationship between 
simulated monthly surface runoff with measured values. In 
general, based on the criteria presented by Moriasi et al., 
(2007), SWAT performed qite well in simulating surface 
runoff for the main stations. 
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Fig. 2 Comparison of the observed (Black line), best 
simulation (Red Line) and 95% prediction uncertainty band 
for surface runoff in different stations. 


Table 1 Monthly model calibration and validation statistics 
for stream discharge 


Station 


Abloo 

Polezokha 

1 

Kerdekhil 

Kiakola 

Righchesh 

meh 

Shirgah 


Calibration Validation 


p-fact 

or 

r-fact 

or 

R 2 

NSE 

p-factor 

r-factor 

R 2 

NSE 

0.12 

0.37 

0.56 

0.55 

0.27 

0.41 

0.55 

0.54 

0.2 

0.34 

0.64 

0.56 

0.38 

0.48 

0.56 

0.48 

0.1 

0.33 

0.44 

0.42 

0.25 

0.37 

0.45 

0.41 

0.06 

0.34 

0.44 

0.43 

0.25 

0.25 

0.41 

0.46 

0.21 

0.59 

0.48 

0.48 

0.36 

0.57 

46 

0.45 

0.16 

0.52 

0.44 

0.43 

0.35 

0.39 

0.48 

47 


Karesang 

0.23 

0.57 

0.58 

0.49 

0.33 

0.56 

0.6 

0.5 

Razan 

0.17 

0.36 

0.62 

0.57 

0.32 

0.43 

0.67 

0.56 
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3.3 Impact of Climate Change on Temperature and 
Precipitation 

Mean annual rainfall for all climate stations during the 
baseline 40-year period (1970-2010) was 731 mm. The 
minimum and maximum Mean annual rainfall in 
Mazandaran province were 578 and 1,307 mm in the east 
and west province, respectively. The average minimum and 
maximum daily temperature were 7/8 and 27.1 C, 
respectively. In Figure 3 the predicted long-term average 
precipitations are compared with the historical data for 
different scenarios. As shown, major changes occur at the 
end of winter and spring, in Mars to June. 




i-T Feb Mar Apr May Jure JiJy A'_g Sep Oct No Dec 


Month 



J-r Feb Mar Ap- May Jure J'Jy A'.£ Sep Oct No Dec 

2 Id nth 

Fig. 3 Comparison of average observed monthly precipitation 
for three GCMs for A2, B1 and A1F1 scenario 

Figure 4 show average monthly changes in maximum and 
minimum temperature for three GCMs and for A1F1, A2 
and B1 scenarios, respectively. T max Increases in temperature 
for A1F1, A2 and B1 scenarios are 2.2, 2.1 and 2.1 °C and 
for Tmin are 2.1, 3.5 and 2.1 °C, respectively. Monthly 
variation in temperature in figure 4 show that maximum 
increase for T max predicted in June and August and minimum 
increase predicted in June and September. Whereas 
maximum increase for T^ predicted in August and 
minimum change predicted in November. In general, all 
projections show an increase in temperature over the basin. 
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Month 


Fig. 4 Comparison of maximum temperatures (left) and 
minimum temperatures (right) for three GCMs and for A1F1, 
A2, and B1 scenarios 

3.4 Impact of Climate Change on Surface Runoff 

Simulation results project a decrease in the annual surface 
runoff of 14.2% in the A2 scenario of CSIRO to an increase of 
21.8% in B1 scenarios of HadCM3 for 2040- 2069. But in 


general, climate change impacts show an increase in surface 
runoff that has a different temporal pattern depending on the 
particular scenario and model (Table 2). The average change 
in annual surface runoff in the main outlets is -1.3%, 5% and 
-1.2% for A1F1, Bl, and A2 scenario, respectively. The study 
conducted by Abbaspour et al. (2009) also reported that 
climate change my increase more frequent and 
larger-intensity floods in the wet regions of northern Iran. 
The monthly variation shows the increase in discharge is more 
pronounced in March and April and the decrease is more 
pronounced at the mid of spring to late summer 
(July-September) (Table 2). In other words, although study 
results show an increase in annual surface runoff, but it 
doesn't occur in a dry season. The increase of surface runoff in 
wet season and decrease in dry season was also concluded by 
Rahman et al. (2012),Yu and Wang (2009), Phan et al. (2011) 
and Shrestha et al. (2013) in different regions. Chang and 
Jung (2010) and Wu et al (2012) also reported that runoff and 
Water yield would increase in spring and substantially 
decrease in summer, respectively. 


Table 2 Predicted relative changes (percent of baseline levels) in monthly surface runoff by different GCMs 


Model - Scenario 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

Ave. 

CGCM2-A1F1 

-2.4 

-1 

47.4 

55.6 

20.1 

17.7 

-3.5 

-13.9 

-12.8 

-9 

-1 

-4.8 

7.7 

HadCM3-AlFl 

-9.9 

16.2 

6.6 

-2.9 

-25.5 

-33.5 

-32.5 

-28.8 

-9.2 

1.1 

-12.1 

-34.3 

-13.7 

CSIR02-A1F1 

3.1 

-4.7 

10.3 

24.7 

1.8 

1.6 

-7 

-14.4 

-17.6 

-4.4 

20.1 

18 

2.6 

CGCM2-B1 

-15 

-3.5 

20.9 

37.1 

21.1 

23.2 

8.8 

0.5 

1.4 

10.8 

0.2 

-14.7 

7.5 

HadCM3-Bl 

8.1 

7 

24 

38.2 

21.7 

25.9 

31 

11.3 

-3.1 

5.2 

16 

11.1 

16.3 

CSIR02-B1 

-1.5 

-6 

3 

2.5 

-20 

-20.7 

-30.2 

-26.6 

-28.6 

-9.6 

20.4 

10.8 

-8.8 

CGCM2-A2 

-3.2 

-1.6 

38.2 

51.1 

18 

16.5 

-2.6 

-14.5 

-11.8 

-6.4 

-0.3 

-8.7 

6.2 

HadCM3-A2 

21.4 

3.4 

12.7 

4.6 

-19.5 

-24.8 

-27.3 

-23.2 

-11.2 

23.4 

60.7 

42.1 

5.2 

CSIR02-A2 

-10.3 

-11.5 

-3.8 

-4.6 

-25 

-27.6 

-32 

-25.7 

-23.8 

-10.5 

-2.5 

-3.4 

-15 


Figure 5 presents the surface runoff Probability for baseline 
and GHGs scenarios for 10 surface runoff classes. Study 
results indicated that climate change may increase high values 
for discharge. As shown in figure 5, the probability of 
occurrence of high values in 10th class (more than 27.4 m 3 /s) 
from 0.3% for baseline has reached to 3%, 2.3% and 2.7% for 
A1F1, A2 and Bl scenarios for a period of 2040-2069. 

Whereas, the probability of occurrence for the most minimum 
surface runoff will decrease in this period. These results 
clearly indicated that climate change will treat water security 
by more floods and severe scarcity and droughts. The increase 
of high values of surface runoff also reported by (Perazzoli et 
al. 2012) in Brazil. 

Fig. 5 Surface runoff probability for baseline and GHGs 
scenarios in different GCMs. 
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IV. Conclusion 

This study assesses the impact of climate change on surface 
runoff in the Mazandaran province basins in the north of Iran. 
To study the effects of climatic variations, the Soil and Water 
Assessment Tool (SWAT) model was implemented to 
simulate the present and future changes in surface runoff. The 
SUFI-2 algorithm in the SWAT-CUP program was used for 
parameter optimization. The climate change scenarios were 
constructed using outcomes of three General Circulation 
Models (CGCM2, HadCM3, and SCIR02) for three emission 
scenarios (A1F1, A2 and Bl). Calibration, validation and 
uncertainty analyses for discharge suggest that the SWAT 
model can be applied to simulate future changes in discharge 
due to climate change. Results indicated that differences 
between the climate models projections in surface runoff are 
high. The study results for 2040-2069 Compared with the 
present climate show an increase and decrease in an annual 
surface runoff with -1.3%, 5% and -1.2% for A1F1, A2, and 
Bl scenarios, respectively. Monthly variation shows that the 
increase in discharge is more pronounced in the wet season 
and the decrease at summer (July-September). The results of 
this study may be helpful to decision makers and other 
stakeholders for adaptive water resource management in a 
changing climate 
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